DNA-condensation, redissolution and mesocrystals induced by tetravalent counterions 
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The distance-resolved effective interaction potential between two parallel DNA molecules is calcu- 
lated by computer simulations with explicit tetravalent counterions and monovalent salt. Adding 
counterions first yields an attractive minimum in the potential at short distances which then disap- 
pears in favor of a shallower minimum at larger separations. The resulting phase diagram includes 
a DNA-condensation and redissolution transition and a stable mesocrystal with an intermediate 
lattice constant for high counterion concentration. 



PACS: 87.15.Kg, 61.20Ja, 82.70.Dd, 87.10+e 

Multivalent polyandries such as trivalent spermidine 
(Spd) and tetravalent spermine (Spe) are abundant in liv- 
ing cells and play a key role in maintaining cellular DNA 
in a compact state [1-3] . They modulate ion channel ac- 
tivities of cells and are essential for normal cell growth. 
Polyamines also facilitate the packaging of DNA in cer- 
tain viruses. The implications can effectively be applied 
in gene delivery and in the field of genetic therapy. Under 
physiologic ionic and pH conditions, the polyamines are 
positively charged and hence DNA is their prime target 
of interaction. 




FIG. 1. (Color online) Typical snapshot in the simulation 
cell. The DNA molecules are shown as two parallel in z direc- 
tion rods over-wrapped by two strings of light grey (neutral 
sphere in MAM (see text), colored yellow in online figure) and 
grey (phosphate sphere in MAM, red in online figure) spheres. 
The tetravalent Spe ions are shown as big black (blue in online 
figure) spheres. Light grey (yellow in online figure) spheres 
represent coions, and dark grey (green in online figure) spheres 
are monovalent counterions. 

In the last decade different experiments have shown a 
condensation and a subsequent redissolution of DNA for 



increasing polyamine concentration C [3-9]. The con- 
densation and redissolution occur at concentrations C c 
and Cd respectively. Between these two thresholds, for 
C c < C < Cd, there is a coexistence of a liquid- like dense 
DNA phase and a very dilute DNA solution. Several the- 
oretical explanations have been presented for the conden- 
sation based on counterion-induced attractions between 
the DNA molecules [1,3,10-15]. However, the origin of 
the experimentally found redissolution transition [10,16] 
is not understood on a molecular level. Proposed mech- 
anisms range from an increased DNA-hydrophilicity in- 
duced by polyamine bindings [5] to Bjerrum pair conden- 
sation [6,7,17] formed by multivalent counter- and mono- 
valent coions and DNA-overcharging [11]. 

In this letter we investigate the condensation and redis- 
solution of DNA on a molecular level by using primitive- 
model computer simulations with explicit tetravalent 
counterions and monovalent salt ions. We trace back 
the condensation and redissolution to the attractions in 
the distance-dependent effective potential U(R) between 
two parallel DNA molecules with R denoting the radial 
distance between their two centers. In fact, the depth 
and position of the attractive minimum play a crucial 
role in whether there is liquid-gas-like phase separation. 
For monovalent microions, the interaction is repulsive 
[18]. For increasing tetravalent counterion concentration 
C, an attractive minimum in U(R) at small separations 
R w 28A shows up which then disappears in favor of 
a shallower minimum at larger separations R ~ 39A. 
Using two-dimensional liquid-state theory for the fluid 
and lattice sums for the solid phases, we calculate the 
phase diagram for columnar DNA assemblies. It includes 
the DNA-condensation and redissolution transition and 
the associated threshold concentrations C c and Cd are 
in agreement with the experimental data. For high con- 
centrations C > 160mM, we predict a stable hexagonal 
mesocrystal with an intermediate lattice constant which 
can coexist either with a dense hexagonal crystal or a 
dilute solution. 

In our computer simulations, we consider B-DNA 
molecules which form a double helix with a pitch length 
of P=34A and Np=20 phosphate charges per pitch us- 
ing the realistic groove-geometry and charge pattern of 
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the Montoro-Abascal model (MAM) [19,20]. A single 
DNA molecule or a pair of parallel DNA molecules, which 
are oriented in the z direction, are placed on the xy di- 
agonal of a cubic simulation box of length L — 102A 
which is three times the pitch length. The box also con- 
tains Nq tetravalent ions, N- — N s + 4Nq monovalent 
coions and N+ = N s + N p monovalent counterions [21]. 
Here N s is the number of added salt ion pairs, and N p is 
fixed by the DNA charge due to the constraint of global 
charge neutrality; N p = 3iVp = 60 for a single DNA and 
N p = 6N P = 120 for two DNA molecules in the box. All 
ions are modelled as charged hard spheres with dq = 8A 
denoting the diameter of the tetravalent counterions. All 
diameters of the monovalent microions are assumed to 
be equal and we choose them to be d c =4A. Hereafter we 
shall call the tetravalent counterions spermine (Spe) since 
experimental data support the idea that it is the charge 
of a counterion, rather than its structural specificities, 
which is important in DNA condensation and reconden- 
sation processes. 

Periodic boundary conditions in all three directions are 
applied. The whole system is held at room temperature 
T = 298K and the water is modelled as a continuous 
dielectric medium with e = 80. The interaction poten- 
tials between the different particle species are a combi- 
nation of hard core and Coulomb potentials. We have 
performed extensive grand canonical molecular dynam- 
ics (GCMD) simulations, as described in Ref. [22], for a 
range of different Spe and salt concentrations. A typical 
configurational snapshot in the simulation box is shown 
in Figure 1. 

First we consider a single DNA molecule in the sim- 
ulation box. It is known that, in the presence of mul- 
tivalent ions, the ionic cloud may not only compensate 
the polyion charge, but even exceed it, resulting in an 
opposite values of the electrostatic potential at some dis- 
tances [17]. The charge compensation parameter of DNA 
phosphate charges, defined as 



ef {-ip Q {r>) + p + {r>)- P -{r>)) 
x 2nr'dr' + N P , 



(1) 



accounts for the integrated total charge at distance r 
away from the DNA axis. Here pi{r) (i — Q, +, — ) are 
radial ion charge densities per pitch length and e > is 
the elementary charge. Data for 9, often called a distance 
dependent effective DNA charge, are plotted in Figure 2 
for different Spe and salt densities. 

For dense Spe concentrations a layered structure ap- 
pears in due to the bulk charge-density oscillations in 
strong electrolytes [23]. The simulations reveal quali- 
tatively different types of competition between the Spe 
and salt ions on the DNA surface depending on whether 
the DNA is overcharged or not. We find that, at a 
constant salt density condition and for C <1.8mM (un- 
dercharged DNA case), the added Spe ions replace the 
small counterions on the DNA surface. For Spe concen- 
tration C >1.8mM (overcharged DNA) i) the influence 



of salt on the parameter 8 diminishes, in other words, 
once the DNA is overcharged, it reluctantly responds to 
added salt, ii) the number of condensed counterions re- 
mains constant, the total adsorbed ionic charge in DNA 
grooves remains constant. It is only the strand ionic 
charge which increases gradually as more Spe is added 
to solution. Thus, the DNA overcharging emerges mainly 
due to the excess Spe charge adsorbed on the DNA phos- 
phate strands. 
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FIG. 2. (Color online) The DNA overcharging parame- 
ter 8 versus the distance r from the DNA axis for different 
Spe concentrations. Thick lines- c s = 25mM, thin lines- c s 
= 100mM. The Spe concentrations are shown next to corre- 
sponding curves. Vertical lines, serving as guide for eye, sep- 
arate the charge layers of different sign for C=400mM. 

The rest of this letter is devoted to the question, how 
the overcharged DNA and the alternating charge layers 
affect the interaction between two DNA molecules. To 
reveal this impact, we calculate the total effective pair 
interaction potential U (R) per unit length for a given 
bulk salt concentration c s and different Spe concentra- 
tions C. The quantity U (R) is gained by integrating 
the distance-resolved interaction force averaged over all 
microion configurations [24]. Results are shown in Fig- 
ure 3. It is seen that even a small trace of spermine 
ions - well below the overcharging threshold - induces an 
attraction between the DNA molecules, except at very 
close distances, see curve (2) for C=0.1mM. This attrac- 
tion has mainly a pure electrostatic origin and arises due 
to charge correlations in the electrolyte. For increasing 
C, this (first) minimum is getting deeper and is achiev- 
ing a maximal depth at the overcharging concentration 
C ~ 1.8mM. There the minimum is mainly resulting from 
entropic forces, i.e. from layering in the Spe number den- 
sity. Further increasing C, again reduces the depth of 
the first minimum. The position of the first minimum, 
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on the other hand, hardly depends on C . The combined 
Spe-layering around the pair of DNA molecules induces 
a second minimum at larger separations as revealed in 
the enlarging inset of Figure 3. This minimum is of elec- 
trostatic origin and occurs for C > 65mM. Again the 
depth of the second minimum increases and decreases 
with C. At intermediate Spe concentrations, we are thus 
confronted with a double minimum potential which is in- 
duced by layering. Bearing in mind that the potentials in 
Figure 3 are scaled for one DNA pitch length, very long 
DNA molecules may easily collapse into the minimum of 
curve (2). This implies that DNA aggregation can take 
place well below overcharging Spe concentrations. 
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FIG. 3. (Color online) Effective pair potential between 
DNA molecules for c a =25mM and C=0mM (1), O.lmM (2), 
0.8mM (3), 18mM (4), 65mM (5), 160mM (6), 280mM (7), 
400mM (8). 

The characteristic double- minimum structure of the in- 
teraction potential U(R) will give rise to unusual phase 
behavior. We have calculated the phase diagram of a 
columnar DNA assembly on the basis of our simulated 
effective pair interactions. We assume that the DNA 
molecules are parallel along a certain length I. This 
length is an additional parameter which we fix to be 
i = 20 x P. We comment on the dependence of the 
phase diagram on I later. The assembly of parallel DNA 
can be considered as a two-dimensional many-body sys- 
tem interacting via U(R) x t and being characterized by 
a DNA particle number density p. We calculated the free 
energies of the fluid and solid phases by using different 
techniques outlined below and perform the traditional 
Maxwell double tangent to identify the coexisting densi- 
ties. 

The free energy of dilute fluid phase is approximated 
by the two-dimensional perturbation theory according to 
Weeks, Chandler and Anderson [25]. The total poten- 



tial is split into a repulsive part U r (R) and a an at- 
tractive part U a (R). The former is identical to U(R) 
but truncated and shifted towards zero at the first min- 
imum at R = R m in- This repulsive potential is then 
mapped onto that of effective hard disks of diameter 
cr e ff using the Barker-Henderson formula [26] o~ e ff = 

a + %n (1 — cxp (— U £^p ))dR. Here the cross-section 

diameter for the DNA molecule is a = 20A. The total 
Helmholtz free energy involves that of a hard disk fluid 

with effective area fraction rj = — and a mean-field 

correction which we simply model as ixp 2 RdR. 
For the former quantity analytical expressions are avail- 
able [27]. The free energy of the solid phase, on the 
other hand, is calculated by a lattice sum assuming a 
two-dimensional triangular lattice as a possible candi- 
date structure. The lattice constant is directly related to 
the DNA number density p. 
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FIG. 4. (Color online) Coexisting DNA densities at dif- 
ferent Spe concentrations C and for c s =25mM. The stable 
phases found are fluid, liquid and two triangular crystals with 
different lattices constants (solid I and solid II). All phase 
transition between these phases are of first order. For the 
sake of better resolution at smaller DNA densities, the ?/-axis 
is expanded below pa 2 < 0.015. 

The resulting phase diagram with the coexisting DNA 
densities is shown in Figure 4 for the whole range of 
Spe concentrations C. At low C there is a strong first- 
order fluid-crystal phase transition. Increasing C widens 
the coexistence region considerably due to the increas- 
ing attractions. Above a threshold concentration of 
C « O.lmM there is enough attraction to stabilize a liq- 
uid phase of high DNA density. The stability of the liquid 
ceases at C w 65mM. At even higher Spe concentrations 
a second crystal with a considerable smaller lattice con- 
stant than that of the high-density solid emerges. We call 
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this novel phase a mesocrystal since its density is inter- 
mediate between that of the fluid and the other almost 
closed-packed solid. 

Another implication of the phase diagram is the con- 
densation and subsequent redissolution. Let us start at 
small C with a dilute DNA solution (see the cross in 
Figure 4) and increase C (dot-dashed line in Figure 4). 
We keep p fixed to pa 2 = 0.002 which corresponds to 
a typical DNA concentration of lmg/ml DNA. First the 
fluid-liquid coexistence line is hit, which implies that the 
system will split into a low density fluid and a high den- 
sity liquid which is the condensation transition. At much 
higher C the coexistence line is touched again and the 
system redissolutes back into the dilute fluid phase. The 
corresponding threshold concentrations of the condensa- 
tion and redissolution are in the range C c « 0.3mM and 
Cd ~ 165mM and agree well with the experimental find- 
ings [10,4]. 

Let us finally comment on the dependence of the phase 
diagram on the DNA length I. Since t is a prefactor of 
the effective potential it plays formaly the role of an in- 
verse system temperature. We have explored the phase 
behaviour for smaller DNA segment lengths of £ = 5P 
and for larger £ = 100P. As a result, the stability of 
liquid pocket depends sensibly on I: it disappers com- 
pletely for small £ but extends towards larger C for larger 
I. The second feature concerns the fluid coexistence den- 
sity: it shifts to considerably higher values for smaller £. 
Hence, condensation and redissolution is prohibited for 
small DNA-segment lengths This is in line with the ex- 
periments of Ref. [28] where a threshold value of £ ss 15P 
for the minimal length £ required for condensation is re- 
ported. 

In conclusion, we have calculated the influence of 
tetravalcnt counterions on the effective interaction and 
the phase diagram of columnar DNA assemblies by 
primitivc-model-type computer simulations and statis- 
tical theories. We find that a small concentration of 
tetravalent counterions induces DNA condensation. The 
layering of the strongly coupled tetravalent counterions 
to the DNA strands yields an oscillatory effective inter- 
action potential with a double minimum structure at 
intermediate counterion concentrations. This explains 
the redissolution transition and triggers a novel stable 
mesosolid. Our threshold concentrations arc in good 
agreement with experimental data. 

Future work should address the phase behaviour of 
DNA solutions which arc polydisperse in their length on 
the basis of the effective interaction found in this paper. 
Furthermore it would be interesting to explore the sta- 
bility of a cholesteric phase [3-7] based on an effective 
interaction which incorporates the relative orientation of 
the two DNA molecules. 
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